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We report our progress in simulating Neuberger valence fermions on Nf = 2 Wilson 0(a)- 
improved sea quarks. We compute correlators with valence quark masses both in the p- and 
in the e-regime, and we match the results with the predictions of the Chiral Effective Theory in 
the mixed regime. This allows us to extract the Low Energy Couplings (LECs) of the Nf ^ 2 
theory and to test the validity of the approach. 

On 
O 



7^. CERN-PH-TH/20 11-229 

> IFT-UAM/CSIC-11-68 

Edinburgh 2011/28 
FTUAM- 11-55 
DESY 11-171 
IFIC/11-54 



The XXIX International Symposium on Lattice Field Theory - Lattice 2011 
July 10-16, 2011 

Squaw Valley, Lake Tahoe, California 



*N.G. acknowledges supported from STFC grant ST/G000522/1 and EU grant 238353 (STRONGnet). 
' Speaker. 



© Copyright owned by the author{s) under the terms of the Creative Coinmons Attribution-NonCoinmercial-ShureAHke Licence. 



http: //pos. sissa. it/ 



Light quark correlators in a mixed-action setup 



Silvia Necco 



1. Introduction 

Lattice QCD simulations with different discretizations in the valence and in the sea sector 
have been performed by several groups and with different setups in the last years [1, 2, 3, 4, 5, 
6, 7]. This is a viable approach to preserve exact chiral symmetry for the valence quarks with a 
"moderate" numerical cost, which could be particularly advantageous for the computation of weak 
matrix elements. 

In our work we apply the Neuberger operator on Nf = 2 non-perturbatively O (a) -improved 
Wilson fermions configurations generated within the CLS common effort. ' In a recent publication 
[8] we have reported our results on the topological susceptibility and the low-lying eigenvalues of 
the massless Dirac operator, which we used to test the so-called mixed Chiral Effective Theory 
approach and finally to extract the chiral condensate. We now extend our analysis to two- and 
three-point functions. 



2. Left current correlator in Chiral Effective Theory 

We consider partially quenched QCD with Ns sea flavors with mass and Ny valence flavors 
with mass m,,, and define the two-point correlators of the left-handed current 

Tr[r'T^]c{xo) = j d\{rQ{x)4{o)), r' = x}fr'Y^p^xif, (2.1) 

where T" is a traceless generator of SU (Ny) and P_ = (1 — Js)/! is the left projector. Low energy 
properties can be described by the con^esponding SU {Ns +Nv) Partially Quenched Chiral Pertur- 
bation Theory (PQChPT); in particular, in a finite volume V = L^T we consider the following two 
situations [9]: 

• niy/s^V » 1: all quarks are in the p-regime. In this case the correlator C{xo) at NLO in 
PQChPT is given by 

, _ fe'ff M„,;effCOSh[M,,;eff(r/2-Xo)] N, dCHO^M^ 

^^""^'-^ sinh[M..,;effr/2] y df ' ^^-^^ 

with 

^eff = ^'(1+Af) = F^jl - ^g'{0;MI) + ^{NsMIu+mIL5)'^ , (2.3) 

<;eff = <(1+Am) = Mif^l + j^[G\0;Mi) + {Mi-Mi)G\0;Mi)] 

- ^ [NsMUu - 2U) +Mi{L5 - 2U)] I , (2.4) 

with the LO pseudoscalar masses M,y = (m, + mj)B and the pseudoscalar propagators in 
finite volume 



'https://twiki.cem.ch/twiki/bin/view/CLSAVebHome 
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j3 = 5.3, csw = 


= 1.90952, 


y = 24^ X 48 




lattice 


K 




amy (e-reg.) 


amv ip-XQg.) 


D4 


0.13620 


0.1695(14) 


0.002 


0.012,0.02,0.028 


D5 


0.13625 


0.1499(15) 


0.002 


0.012,0.02 


D6 


0.13635 


0.1183(37) 


0.002 


0.006,0.015,0.025 



Table 1: Simulation parameters 



• myEV < 1, mjy » 1: valence quarks are in the e-regime, sea quarks are in the p-regime. 
In this case the NLO left correlator at fixed topology is given by 



with /ly = m^iy and 



IT 



2T^ 



NsdG\Q;Ml/2) 
2 dT 



/II (T) 



F^(l+Af|M„,=o 

2 



^^G\0-Mi/2) + ^{NMlU: 



1 



Ov{x) = x[Iv{x)Kv{x) +/v+i(x)A'v-i(x)] + 



(2.6) 



(2.7) 

(2.8) 
(2.9) 



where ly, Ky are modified Bessel functions. 

A well-established strategy to extract the LO {F,B) and the NLO (L,) Low Energy Couplings 
(LECs) is to match lattice results for the two-point correlators with the predictions of the ChPT. 
Depending on the chiral regime, the systematic uncertainty associated to the matching can vary 
considerably. Exploring different regimes in a finite box can then be useful to give reliable esti- 
mates of systematic errors associated to the LECs. 



3. Lattice results 



3.1 Setup and results for the Low Energy Couplings 

We have carried out our simulations on CLS lattices of size V = 24^ X 48, with jS = 5.3 
corresponding to a lattice spacing a = 0.0649(10) fm [10]. We have considered three p-regime 
sea quark masses k = 0.13620,0.13625,0.13635, labeling the sets by D4,D5,D6 respectively. On 
these configurations we have built the Neuberger-Dirac operator, with valence quark masses both 
in the p- and in the £-regime. The simulations parameters are collected in Table 1, where we report 
also the effective sea pion masses in lattice units (see [8] for further details). For valence quark 
masses in the p-regime, we have collected ~ 80 — 160 measurements. For the £-regime, we have 
in total ~ 70 — 100 measurements of the left correlator in the topological sectors | v | =0,1,2,3. 

We have computed the two-point left current correlator defined in Eq. (2.1) adopting the low- 
mode averaging technique of [11]. For the renormalization of the local axial cuiTcnt we used 
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Figure 1: Left: the renormalization constant of the local axial current extracted through the chiral Ward 
identity at fixed topological charge saturated with the zero modes. The blue,green and red points refer to the 
lattices D4,D5,D6 respectively. Right: afg eff extracted from the left correlator in the mixed p/e-regime as a 
function of the topological charge. The gray bands represent the average over |v| — 1,2,3. 



the chiral Ward identity at fixed topological charge saturated with the zero modes, following the 
strategy proposed in [12]. In particular we considered 

Zr \im {m,y)Y^{d,P"{x)jm)v = lim (m2y)£(P''WP-(0))v, (3.1) 

X X 

where P" = i//r"75i/A is the pseudoscalar density. The right-hand side is completely saturated with 
zero modes and can be directly computed without chiral extrapolation; on the left-hand side, the 
chiral limit has to be taken numerically. The renormalization factor Z™" can then be extracted 
from ratios of correlations functions. For those correlators we have collected a total of ~ 70 — 140 
measurements for the topological sectors < | V | < 5. On the left side of Fig. 1 we report the results 
for the different sea quark masses corresponding to D4,D5,D6 and for v = 1,2,3,4. Apart from 
some fluctuations in the topological sector |v| = 1, the data do not depend on v and on the sea 
quark masses. From this study we obtained 

Z7'' = 1.65(1). (3.2) 

We then extracted Mvv;eff, /^ff, Fe;eff from the left correlators, according to Eqs. (2.2, 2.6); on the 
right side of Fig. 1 we show for instance the (bare) pseudoscalar decay constant computed from the 
mixed p/£-regime correlators at |v| = 0, 1,2,3. As it should be, results do not depend significantly 
on the topological charge; fe;eff has been determined by averaging over |v| = 1,2,3 (gray band 
in the plot). We then performed a global fit of pseudoscalar masses and decay constants using 
the NLO predictions of the PQChPT given in Eqs. (2.3,2.4,2.7). From this fit we obtained the six 
parameters 

af'^ = 0.0259(10) (8) , aS^ = 1.52(11)(9), 

Lf = 0.00054(10) (4) , Lf = 0.00073(15)(2), (3.3) 
(L4-2L6)'^ = 0.00033(5) (0) , (L5 -2L8)^ = -0.00025(15)(5), 

where "R" means that the couplings are renormalized. The constants Lf are given at the scale jX = 
110 MeV, while is given in the MS scheme at /i = 2 GeV. For the tree-level pseudoscalar masses 
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am am 

V V 

Figure 2: Global chiral NLO fit for M,,,,;eff, ^eff- The squares in the right plot correspond to the e-regime 
results, Fg;eff. 



we used M,y = (mf + mJ)B^; for the sea Wilson quark masses we used the PC AC masses computed 
in [13, 10], with renormalization constants Zj = 0.78(1) [14], Zf(MS, 2 GeV) = 0.543(8) [15]; 
for the overlap valence quarks we have computed Z™''(MS, 2 GeV) = 1.88(15) like in [8]. The 
first error is statistical, while the second one reflects the uncertainty on the PCAC quark masses. 

The global fit is shown in Fig. 2. It is important to stress that having valence quark masses in 
the £-regime allows us to extract the pseudoscalar decay constant directly at amy = and hence to 
have a good control over the chiral extrapolation in m,.. 

Appropriate combinations of the L, yield the conventional SU(2) LECs Ti; in particular we obtain 

Z3 =4.46 (30) (14), [4 = 4.56(10) (4). (3.4) 

The FLAG working group [16] gives an estimate fi, = 3.2(8) and quotes I4 = 4.66(4) (^j'*^) [17] 
as best determination of I4. Our determinations are in agreement with those results within the 
uncertainties. By converting aF^ and aB^ in physical units we obtain = 79(4) MeV and 
r'^(MS, 2 GeV) = B^(F'*)2 = (306(1 1) MeV) 3, which also he in the right ballpark. 

3.2 Unitary point 

A natural definition of the unitary point corresponds to the matching of the pseudoscalar 
masses 

■^vv;effln!,.=m,;un = ^ss;eff- (3-5) 

We have extracted mv;un through a linear interpolation of M^^,.^^^. The pseudoscalar decay constants 
in the full Wilson theory (^/gff) and in the mixed action theory at my-un (^effl'"i ="iv un) should then 
be equal up to discretization effects. In Fig. 3 we show F^\,n,,=m^.„^ and F^^f^ for the three lattices 
D4,D5,D6. The Wilson on Wilson data have been taken from [13]. Within our statistical accuracy, 
we do not observe large cutoff effects close to the unitary point. 

Discretization effects in a mixed action setup can be investigated by means of Mixed Action 
PQChPT [18]. If all mass scales are in the p-regime and one assumes a power counting a ~ v ~ 
{generically small masses), it is easy to see that "unitarity violations" of 0{a^) appearing in M^^ 
and {MI^ — M^^) affect our observables only at NNLO, and hence are expected to be small. Our 
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Figure 3: Effective pseudoscalar decay constant computed in the full Wilson theory (black points) and in 
the mixed action theory at the unitary point (red points) as a function of the pseudoscalar mass. The empty 
square represents the decay constant in the chiral Umit extracted from the chiral fit, Eq. (3.3). 



numerical data seem to confirm this expectation, although a more detailed study is necessary in 
order to have a precise estimate of those effects. ^ 



4. K TTTt amplitudes and AI = \/2 rule 

In [19], a strategy to study the role of the charm quark in the A/ = 1/2 rule has been proposed. 
The basic idea is to start from the GIM limit m„ = nid = nis = nic, write down the AS = 1, CP- 
even effective weak Hamiltonian and match it with the effective Hamiltonian in the SU(4) Chiral 
Effective Theory. At leading order in ChPT, the ratio of amplitudes corresponding to A/ = 1/2 
and AI = 3/2 K ^ nn transitions can be expressed in terms of given LECs gf. As usual, those 
couplings can be computed by matching suitable correlation functions computed on the lattice with 
their analytic expressions obtained in ChPT (see [19] for details). Quenched results obtained with 
Neuberger fermions have been presented in [20]; this strategy, for which exact chiral symmetry 
at finite lattice spacing is an essential ingredient, can be implemented also with our mixed action 
setup. Preliminary results for the bare ratios of three- and two-point functions defined in [20] 
have been obtained so far for the lattice D4 (interpolated to the unitary point) 

/?+= 0.64(6), R- = l.5{2), (4.1) 

of the same order as the quenched results for pseudoscalar masses in the same range. 
The second part of the strategy consists in studying the amplitudes ratio as a function of > 
/M„ = m^i = nis and test if by increasing the charm quark the A/ = 1/2 enhancement observed 
experimentally is reproduced. 



5. Conclusions and outlook 

Mixed action simulations with exact chiral symmetry in the valence sector provide a powerful 
tool to investigate QCD properties. Given the significant errors associated to chiral extrapolations, 

recent study with Neuberger valence fermions and Wilson Twisted Mass sea quarks [6] indicates large cutoff 
effects close to the unitary point. 
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it is important to explore several chiral regimes where the systematic effects are different, so that 
such errors can be properly quantified. In particular, we have considered valence quarks both in 
the p- and in the £ -regime and matched the lattice results with the PQChPT predictions at NLO. 
A global chiral fit of pseudoscalar masses and decay constants allowed us to determine F , B, I3, 14 
with comparable accuracy with respect to other results present in the literature. By matching pseu- 
doscalar masses in the valence and in the sea sector, we did not observed large cutoff effects close to 
the unitary point. We stress that a precise determination of the LECs would require an improvement 
of the statistical precision and a careful monitoring of autocorrelations; furthermore, a larger set 
of volumes, lattice spacings and quark masses would guarantee a solid estimate of the systematic 
uncertainties associated to lattice artefacts, finite- volume effects and chiral extrapolation. 

As next step, we are implementing the strategy proposed in [19] to investigate the role of 
the charm quark in the A/ = 1/2 rule. This calculation is more challenging and requires further 
computational efforts. 
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